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One of the major differences between fermions and bosons is that fermionic states have a max- 
imum occupation number of one, whereas the occupation number for bosonic states is in principle 
unlimited. For bosons that are made up of fermions, one could ask the question to what extent the 
Pauli principle for the constituent fermions would limit the boson occupation number. Intuitively 
one can expect the maximum occupation number to be proportional to the available volume for 
the bosons divided by the volume occupied by the fermions inside one boson, though a rigorous 
derivation of this result has not been given before. In this letter we show how the maximum oc- 
cupation number can be calculated from the ground-state energy of a fermionic generalized pairing 
problem. A very accurate analytical estimate of this eigenvalue is derived. From that a general ex- 
pression is obtained for the maximum occupation number of a composite boson state, based solely 
on the intrinsic fermionic structure of the bosons. The consequences for Bose-Einstein condensates 
of excitons in semiconductors and ultra cold trapped atoms are discussed. 



Textbooks on quantum mechanics teach that fermionic states have a maximum occupation number of one, whereas 
the occupation number of bosonic states is unlimited. However, if the boson is actually a many-fermion state (e.g. an 
atom), the Pauli principle results in a maximum occupation number for the bosonic state, irrespective of the nature 
of the boson-boson interactions. Though one can easily argue that the maximum occupation number for a composite 
boson state should be of the order V/r D , with V the available volume, D the dimension and r the range of the 
intrinsic fermionic wave function of the boson Q| , a rigorous and accurate derivation of this limit has not been given 
before. Here we try to quantify this effect directly in terms of the internal fermionic structure of the boson. We find 
that the Pauli principle puts an upper limit for the occupation number of a composite boson state, irrespective of the 
boson-boson interactions. For a boson state composed of two fermions, this maximum occupation number (MON) can 
be calculated as the maximal eigenvalue of a fermionic generalized pairing operator. Using variational principles we 
can derive a very accurate analytical estimate of this eigenvalue. This allows us to write down analytical expressions 
for the MON in a number of cases: for excitons in semiconductors, for Gaussian wavefunctions and for atoms in 
Harmonic traps. 

II. THE MAXIMUM OCCUPATION NUMBER AND THE GENERALIZED PAIRING PROBLEM 

As a first case, consider a number of bosons occupying a given one-boson state V'o- The operator that creates one 
boson in the state tpo can be written as combination of pairs of fermion creation operators: = • Bc^jjaffij- The 
number of bosons occupying the state i/>o can be evaluated as the expectation value in the many-boson state of the 
boson number operator Nq = SJ&o- The MON is then given by the maximal eigenvalue of the number operator No. 
According to the Bloch-Messiah theorem we can always find a unitary transformation that brings the matrix Bq in a 
canonical form: 



where the indices (a, a) correspond to conjugated pairs of canonical fermion states. After this transformation, we see 
that the boson number operator is equivalent to a fermionic generalized pairing operator 
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I. INTRODUCTION 
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The MON corresponds to the ground-state energy of the pairing Hamiltonian H = —G. This energy can be calculated 
efficiently using algebraic methods ||||f§, for particle numbers up to a few thousand. Here, instead, we derive a very 
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accurate and practical analytical estimate of the MON. At the same time, this gives a simple and accurate analytical 
estimate of the ground-state energy of the pairing Hamiltonian H . 

Upper bound: Consider the matrix representation of the operator Nq in the basis spanned by the states that are a 
product of N pairs of fermions (a, a). The eigenstate corresponding to the largest eigenvalue is contained in this 
space fa]. The Gershgorin circle theorem results in the following upper bound: 



MON < max 




(3) 



where S = ^ Q \c a \, a denotes a state of N pairs distributed over the canonical pair states (01, . . . , ctjv), and the c a . 
denote the corresponding Bloch Messiah coefficients. Because of the normalization of the fermion pair state of Eq. ([!]), 
the first term in Eq.(||) is bounded by one, 

N 

£M 2 <EM 2 = 1 - (4) 

3=1 a 



The second term is bounded because 



x(S - x) < 5 2 /4, (5) 



for any real number x, in casue for x — Y^j=i \ c ^j I- Using these facts we can obtain a slightly relaxed but more general 
bound, valid for any N: 

MON < — + 1. (6) 
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Lower bound: Taking the following fully paired state as a trial state, 

\MN))= (l>^) 1°)' ( ? ) 
where |0) denotes the zero-particle vacuum state, we can derive a lower bound from the variational principle that 

MON > {ip T \G\ip T )/{^T\^T}- (8) 

This results in 

N\N-1 + (M - N)S 2 ] 

MON> — L - — r — - — i-, (9) 

M(M - 1) w 

with N the number of bosons (fermion pairs) in the trial state and M the number of canonical pair states o^at|0). 
The maximum of the lower bound is found to be very close to the upper bound. For S 3> 1, this maximum occurs at 
N ~ M/2. Taking N — (M + l)/2 in Eq.@ leads to 

MON>i[5 2 + l](l + ^). (10) 

For a realistic description we will have to take the limit of M going to infinity, because the Hilbert space has an 
infinite dimensionality. This means that the total number of fermion pairs can go to infinity too (the quantity N in 
Eq-© )• If the quantity S has a finite value in this limit, then the upper bound in Eq.(|6j) is finite, and hence also the 
lower bound in Eq.(|9[) will have a finite value. Eq.(|^) and Eq.(|lO|) allow one to put general but very stringent limits 
on the MON: 

^- + 7 <MON< ^- + 1. (11) 
4 4 - - 4 

To estimate the MON we have to calculate the quantity S. In two cases we can derive detailed expressions for this 
quantity: in the case that the one-boson wave function corresponds to a plane wave (for instance in a cubic volume 
V with periodic boundary conditions) and in the case of a Gaussian one-boson and intrinsic wave function. 
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III. THE MAXIMUM OCCUPATION NUMBER FOR AN EXCITON STATE 



For a boson in a cubic volume V with periodic boundaries, we can express the wave function in the momentum 
coordinates: 

V(kx,k 2 ) cx 5(ki - k 2 )F [(k x + k 2 )/2] , (12) 

where F(k) corresponds to the Fourier transform of the real space intrinsic wave function. The canonical pair states 
of the Bloch-Messiah decomposition for such a state are given by the momentum pairs (k, — k), 

6S = X>(k)a£ciU (13) 

k 

Then the quantity S from Eq.(|ll]) can be calculated as 

S = £|F(k)|. (14) 

k 

This expression can be applied directly to electron-hole excitons in semiconductors 0J^,^,||- Two recent let- 
ters |l(],[ll| addressed the question to what extent these excitons can be regarded as bosons. A criterion to assess 

that, was derived from the expectation values and \B,Bn \(/> iN) ), with & the one-exciton 

creation operator and 



\^)oc(m |0). (15) 



They expect the boson picture to be valid up to 



Nsi -k> (16) 



with V the available volume, a x the exciton Bohr radius and A a scaling factor. Their model allows them to propose 
a value A ~ 100 up to 400. Compare this with a value of A ~ 1 that was suggested from the disappearance of 
the exciton binding energy at the Mott density, due to Coulomb screening g§]. The state \<j>^) can be seen as 
a number-projected Hartree-Fock-Bogoliubov state. It can be expected that the true many-exciton state will differ 
significantly from this mean-field state at high exciton densities. Whatever the true state might be, the occupation 
number B) will be limited by the upper bound for the MON from Eq.(|^). Assuming that the intrinsic wave 
function of the excitons takes the form of a Is hydrogen wave function, we obtain the momentum representation 



with a x the exciton Bohr radius. We find that S 2 — V/(Tra^.), such that 

MON=^. (18) 

This means that the Pauli principle forbids having more than V/^Tra^.) excitons. The corresponding scaling factor A 
is given by 

A = 4tt. (19) 
For the two-dimensional box with surface V, we have 

(l + Pa2/4) 3/2 

and hence MON = V/(^a 2 x ) and A = tt/2. 
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IV. THE MAXIMUM OCCUPATION NUMBER FOR A GAUSSIAN STATE 



Another case where we can derive an expression for the quantity S, is in the case of Gaussian wave functions for the 
one-boson state and for the intrinsic two-fermion state. For that it is useful to work with a reduced density matrix, 
integrating out one of the fermion coordinates: 



P( r i.r'i) = J ^(ri,r 2 ) t V(ri,r2)dr 2 . 



(21) 



The values c a of the Bloch-Messiah decomposition correspond to the square roots of the eigenvalues of the matrix p. 

Assuming a Gaussian shape for the wave functions we can write the reduced density 



So we can write S 
matrix as 



Tr 



, a <- 2 +' J . 2 _£ 

p (r, r ) oc e w 2 e 



r') 2 
72 



(22) 



where W gives the width of the one-boson state, while the scale v corresponds to the width of the fermion distribution 
inside one boson. Note that the square of a matrix of this form has again the same form, with different coefficients. 
By equating these coefficients, we can find an explicit form of p? , and hence the trace, 



S 



This means that we can estimate the MON by 






3/2 






2^" 












V 



.3/2 



(23) 



(24) 



The cubic exponent in Eq.(|lq) and Eq.(p4[) demonstrates that there is a volume effect: for every condensed boson a 
minimal volume is required that is proportional to a 3 or v 3 . 



V. THE MAXIMUM OCCUPATION NUMBER FOR A TRAPPED ATOM 



For the derivation of Eq.(|24|) we assumed Gaussian wave functions. The wave function of an atom in an harmonic 
trap will closely resemble this form for the center-of-mass coordinate. The intrinsic wave function will not be a 
Gaussian. Still we can derive a useful estimate of the MON using a Gaussian approximation. To set the case we 
consider hydrogen atoms in a harmonic trap. Let us define the electron overlap function f(r — r') as 

f(r r') = J 4nt(r r>i„t(r" - r')dr" 

= (^nt|e (r - r ' ) ' V |^nt), (25) 

where ipint denotes the internal wave function of the atom. The function /(r) is equivalent to the Fourier transform 
of the momentum distribution. If we can assume that the overlap function is nearly Gaussian, i.e. if 

/(r)~e-&, (26) 



then the density matrix will still have the form of Eq.(22), so that Eq.(g4J) applies. Note that t he first order term 
in a series expansion of f(r) vanishes for any differential; 



ible wave function ipint , and therefore Eq. ( |26| ) has the correct 
leading order. For the Is orbital of the hydrogen atom, the overlap function is given by /(r) = /(r/ao), with ao the 
atomic Bohr radius and 

/(*)= (l + s + y) z~ X - (27) 

Using a series expansion we can write 

2 4 5 

■^w— T + H-a + -- < 28) 
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Keeping only the leading term, we can write 



/(r)~e 6 "o. (29) 
This is an approximation, but because the next-to- leading order term is of order r 4 , we believe that it is good enough 



to allow for a reliable order-of-magnitude estimate of the MON. The function .F(r) in Eq.(29) has the Gaussian form 
of Eq.©, with 

v ~ -^a . (30) 

The MON becomes 

16W 3 

M ° N =3VR' (31) 

Another way to arrive at this result, is based on the relation of the displacement operator x ■ V to the kinetic energy 
T, valid for s electrons : 

(lp int \ (r • V) 2 |V>t„t) = -(V'tntl ( r ' P/^f \i>int) 



2r 2 m e T 



(32) 



with m e the reduced mass of the electron. From the virial theorem we know that the kinetic energy T is equal to the 
electron binding energy Ef,. This allows to estimate v as 

^ (33) 



^/8m e Eb 



It is easily checked that Eq.(|30D and Eq. (|33[) coincide for the hydrogen Is orbital, while the latter is a bit more general. 

Spin degrees of freedom can be taken into account straightforwardly: for hydrogen and the alkali atoms, where the 
valence electron has one spin-1/2 degree of freedom, the estimate of Eq. (p!T[) for the MON has to be multiplied by a 
factor of 2. 

For atoms built up of more than two fermions, the number operator is no longer equivalent to a fermion pairing 
operator. However, we can still apply the above formalism by considering one valence electron and the singly ionized 
atom. This makes sense for the alkali atoms, because the wave function of the first electron out of a closed shell has 
a very long tail and a small overlap with the other electrons. In this case, Eq.(^3|) can be used again to obtain an 
order-of-magnitude estimate, with E\, replaced by the first ionization energy. Because of the exchange between valence 
and core electrons, it is no longer guaranteed that this will result in a strict upper bound for the MON. However, it 
is still useful as an order-of-magnitude estimate. To find a better estimate for the MON in these cases will require 
advanced numerical many-body techniques. 

Applying Eq.(|3^) to the singly ionized atom and the valence electron, we can obtain an estimate for the MON from 
the experimental electron binding energies ]l3|]. F or a trap of width l/i we find a MON of the order of 10 13 , for all the 
alkali atoms. Present day experiments p^ , |15| , |16| , p7[ find occupation numbers of the order of 10 6 - 10 9 , which means 
that they are in a regime of very low densities, such that Pauli effects due to the constituent fermions can safely be 
neglected. 

At higher densities one can expect serious alterations in the condensate wavefunction long before the MON is 
reached, due to the interatomic interactions. However, even at high densities the above formalism remains valid, 
provided that one considers the right one-boson state. As presented in his famous 1962 paper |l8[ , C.N. Yang defines 
the wave function of a Bose condensate as the eigenvector with the largest eigenvalue of the one-body density matrix. 
This Bose wave function can be decomposed as in Eq. ([!]). If a condensate would persist even at high densities p9| , 
its occupation number remains bounded by the expression for the MON derived above. As a matter of fact, nuclear 
physics might offer an example of this situation. Recent large-scale shell model calculations |^,^l[ indicate that the 
valence neutrons of open-shell nuclei tend to form bosonic pairs, and that these pairs might have high occupation 
numbers (compared to the number of single-particle states in the system) , that might be very close to their maximum 
value. 
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VI. CONCLUSIONS 



We conclude that the Pauli principle for the constituent fcrmions puts an upper limit to the occupation number 
of a composite boson state. For electron-hole excitons in semiconductors, this maximum occupation number can be 
calculated by solving a fermionic generalized pairing problem. A very accurate approximation to the exact solution 
shows that there is a volume dependence in the maximum occupation number: each exciton consumes a volume that 
is proportional to the extent of the electron- hole intrinsic wave function. We find a maximum of V/(4Tra^,) excitons in 
a three-dimensional volume V ( V/ (f-a^) excitons in two dimensions). The expression also puts an upper limit to the 
maximum occupation number of an atomic Bose Einstein condensed state. A rough estimate based on the electron 
binding energies, leads to values around 10 13 for a trap of width lfi, much higher than the occupation numbers found 
in present day experiments. This means that we can safely look at these condensates as truly bosonic condensates, 
though we know that actually they are made up of fermionic particles (nucleons and electrons). 
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